Longitudinal versus transverse spheroidal vibrational modes of an elastic sphere 
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Analysis of the spheroidal modes of vibration of a free continuum elastic sphere show that they can 
be qualitatively grouped into two categories: primarily longitudinal and primarily transverse. This is 
not a sharp distinction. However, there is a relatively stark contrast between the two kinds of modes. 
Primarily transverse modes have a small divergence and have frequencies that are almost functionally 
independent of the longitudinal speed of sound. Analysis of inelastic light scattering intensity from 
confined acoustic phonons in nanoparticles requires an understanding of this qualitative distinction 
between different spheroidal modes. In addition, some common misconceptions about spheroidal 
modes are corrected. 



PACS numbers: 62.20.-x,43.20.Ks,62.25.+g 

I. INTRODUCTION 

With the explosion of interest in the optical properties 
of nanoparticles, the classic elastic mechanical problem 
of the vibrational modes of a free continuum sphere has 
found a new context for application. The problem was 
formally and numerically solved back in 1882 J» Nanopar- 
ticles, i.e. spherical clusters of atoms ranging in diameter 
from 1 nm to 100 nm, have sufficiently few atoms that 
the continuum approximation can be questioned^^ Even 
so, it is acceptable to ignore the effects of the discrete 
crystal lattice for the few vibrational modes with lowest 
frequency as long as the nanoparticle diameter exceeds 
several nanometers. 

Inelastic light scattering of a continuous laser beam 
shining on the nanoparticle permits detection of the me- 
chanical vibrations since the changing size and shape 
of the nanoparticle modulates the polarizability of the 
nanoparticle, so that the monochromatic incident light 
turns into scattered light with sidebands shifted up and 
down by the frequency of the vibrations. This can be 
seen using experimental setups of the Raman and Bril- 
louin type. 

For theoretical convenience the material is assumed to 
be homogeneous, isotropic and linear. The outer sur- 
face of the sphere is free from externally imposed stresses 
and this situation will be referred to as the "free sphere 
model" (FSM). The original paper by Lamb 1 classified 
the FSM modes of vibration into two classes, now called 
"torsional" (TOR) and "spheroidal" (SPH). The distinc- 
tive feature of torsional modes is that the material den- 
sity does not vary. In other words, the divergence of the 
displacement field is zero. Furthermore, the spherical 
symmetry permits classification of the modes by angular 
momentum number £ > 0. (However, later on we will 
show that there is value in considering £ to be a continu- 
ous variable.) There is no dependence of the frequency on 
the z angular momentum to. Beyond this, the modes are 
indexed by n > 0. It is convenient to let p denote either 



SPH or TOR, to indicate individual modes by (p, £, n) 
and their frequencies by Lu p i n . 

What we explore in this paper is an additional classi- 
fication of the SPH modes beyond that employed since 
Lamb. In particular, SPH modes can be classified (albeit 
approximately and subjectively) as either being primar- 
ily longitudinal (SPHl) or primarily transverse (SPHy) 
in nature. The specific meaning of this will be explained 
further on. This is not a sharp division, and actual modes 
fall somewhere in between the two ideals. However, the 
contrast is sufficiently sharp that this new distinction 
among SPH FSM modes as SPHl or SPHt is a very 
important tool. 

In a recent theoretical paper, G. Bachelier and 
A. Mlayahi predicted that (SPH,£ = 2,n) modes with 
differing values of n contribute to the Raman spectrum 
in a highly non-uniform way. In this paper we will show 
that this can be explained using the previously men- 
tionned distinction between SPH modes. They pointed 
out that there are two separate mechanisms that cou- 
ple (SPH/ = 2) acoustic vibrations to the surface plas- 
mon resonance and in turn lead to Raman scattering. 
First, change of the particle shape and second, modula- 
tion of the density leading to change of optical properties 
through the deformation potential. 

Section[H]reprises the formalism necessary for the FSM 
solution. In Section lTTTl we show explicitly what we mean 
by SPH L and SPH T . In Section |IV| we illustrate the 
natural appearance of SPH^ and SPHt modes in the 
high frequency limit. Section discusses these results 
and their connection with inelastic light scattering ex- 
periments. 



II. THE FREE SPHERE MODEL 

Vibrational modes of a free linear elastic continuum 
homogeneous isotropic sphere were found by Lamb in 
1882. 1 
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For a mode with angular frequency u, the displace- 
ment of material point f from its equilibrium position 
is u(r)cos(wt). For a m = TOR mode, u = A 
V x (rji(kTr)Pg{cos9)), where ji are spherical Bessel 
functions of the first kind and Pg are Legendre polyno- 
mials. For a m = SPH mode, u = ul + ut where 



u L (r,0) = B\7j e (k L r)P e {cos9) 



(1) 



and 



ur(r, 9) = CVxVx (fj e (k T r)P e {cos9)) (2) 

where A, B and C are real coefficients, vl^l = vtIzt 
= to, and i;t and vl are the transverse and longitudinal 
speeds of sound. 

Modes with z angular momentum m ^ have a differ- 
ent functional form. 

R is the nanoparticle radius. If cr^ is the stress tensor, 
the boundary conditions at r — R are a rr = a r g = 0. It 
is convenient to introduce dimensionless frequencies r\ = 
krR and £ = UlR- Following Eringen, 5 application of 
these boundary conditions determines the allowed SPH 
vibrational frequencies as zeroes of a 2 x 2 determinant 
for I > 0. 
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For £ = 0, the allowed vibrational frequencies are the 
zeroes of Tn. 

Noting that the displacement fields are real-valued, it 
is appropriate to use the following inner product between 
two displacement fields ua and ub& 



(u A \u B ) 



J r<R u A {r) ■u B (/)pd 3 f 

Ir<RP d3ff 



(4) 



A normalization condition (such as (u\u) — 1) would 
typically determine the final values of B and C . But 
the details of the condition do not affect the results re- 
ported here. The displacement field u{r) for some se- 
lected modes are depicted in Fig. ^ 



III. SPHEROIDAL MODE LONGITUDINALITY 

Isotropic elastic materials differ in their Poisson ra- 
tio, v, which is related to x — vt/vl through x = 
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(SPH,2,2) 
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FIG. 1: (Color online) Displacement fields u(r) of selected 
SPH 1=1 modes. As explained in the text, the first three 
are primarly transverse (i.e. SPHt). (SPH, 2, 3) is primarly 
longitudinal (i.e. SPHl). The equilibrium surface of the 
nanoparticle and the z-axis are shown as dotted lines. The 
solid (red online) line shows the distorted surface. Note that 
the (SPH, 2,1) mode does not change the nanoparticle shape. 



y(T— 2v)/(2 — 2v). Figure shows how the dimension- 
less frequency, rj, of the SPH I — 2 FSM modes varies 
with vt/vl- It is quite apparent that some modes keep 
the same r\ as vt/vl is varied. However, other modes 
change frequency as vt/vl changes. There are transition 
points where a given mode changes from being constant 
to varying with vt/vl- 

This pattern visible in Fig.[5Jmotivates the search for a 
numerical criterion to permit this contrast among modes 
to be quantified. We adopt the starting point that in 
some sense some modes are more transverse in nature 
(SPHt) while others are more longitudinal (SPH^). We 
then coin the term "longitudinality" , denoted by L, for a 
quantity that varies on a scale from to 1 with being 
purely transverse and 1 being purely longitudinal. There 
is no single obvious way of doing this. Rather, we have 
evaluated a number of quantities as candidates for the 
best measure of longitudinality, of which we present four 
which work well. These will be denoted LI, L2, L3, and 
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FIG. 2: (Color online) Dimensionless mode frequency r\ as a 
function of vt/vl for (SPH/ = 2) modes. Vertical lines (blue 
online) mark vt/vl for Au, Ag, Si and Ge from left to right. 



LA. 

Consider a particular SPH mode with indices 
Its frequency is lu(vl,vt). Define LI by 



and 



LI = 



vl duj 
w dvL 



x dr/ ^ 
rj dx 



Vt duj 



(5) 



Noting that u(r) = urir) + UL(r), we define LI as 
(ul\ul) / (u\u) , and L3 as 1 — ((ut\ut) / (u\u)) . But note 
also that (ul\ul) + (ut\ut) ^ (u\u) since (ul\ut) i= 0. 

Given a fixed value of vt /vl and n, rj may be consid- 
ered to be a continuous function of t, as in Fig. [3] In 
terms of this rj(£), define 



L4 = VT 

vl — vt \ir d£ 



(6) 



Let < .. >y and < .. >s denote averages over the 
nanoparticle volume and surface, respectively. In partic- 
ular, (u\u) =< Ur + Ug + u^ >v- Some other measures of 
interest are as follows: URV = < v? r >y /(u\u), URS = 
< u 2 r > s /(u\u), UTS = < u 2 e + u\ > s /(u\u), and U2S 
= URS + UTS. 

Note that all of these quantities are defined in such a 
way as to be independent of m. 

Except at low r/, Fig. 0] shows that LI and L2 are in 
close agreement. L3 and L4 are not plotted, but also 
agree closely except at low rj (see Table . Generally, 
a given mode either has all of LI, L2, L3, and L4 low, 
or else all high. It is thus possible to classify modes as 
SPHi or SPHt- Rarely, there are cases where the values 
of LI, L2, L3, and L4 are in the intermediate range, such 
as in Fig. [3] for £ = 3 for the two modes near r\ = 13. 
Such modes which are neither clearly SPH/, nor SPHr 
always occur in pairs. The reason for this is explained in 
Section El 

Table [H] provides additional information about the 
modes. The dimensionless frequency r/ is provided for 
convenience, as is the ratio of coefficients B and C. 



FIG. 3: (Color online) Variation of the SPH mode frequency 
with i for a material with vt/vl = 0.5. Full circles are exact 
FSM frequencies and are connected with curves calculated for 
non-integer I (red online). Lines with crosses are roots of Tn 
(blue online) which approximate SPH_l modes. Lines with 
empty squares are roots of T43 (green online) which approxi- 
mate SPHt modes. 




FIG. 4: (Color online) Dimensionless mode frequency r\ as a 
function of longitudinality measures LI (circles (red online)) 
and L2 (crosses (green online)) for SPH modes of a material 
with vt/vl = 0.5. t = 1 . . .4 from left to right. SPHt modes 
have LI and L2 near zero, while for SPHl modes they are 
close to 1. 



In principle, C/B could be expected to provide useful 
information about whether a mode is SPED or SPHt- 
In the extreme case that C = 0, the mode is evidently 
SPHl , and likewise when B = the mode is SPHt ■ But 
the C/B values do not exhibit an informative pattern. 

There is a strong contrast in the values of URS for dif- 
ferent modes. However, it does not correlate to whether 
the mode is SPH^ or SPHt except at high enough rj. 
URS is an interesting quantity because it is the one we 
have to monitor for the surface deformation mechanism 
except for £ = modes. 
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TABLE I: Longitudinality measures for a material with vt/vl 
= 0.5. 
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Group theoretical arguments^ show that only SPH 
modes with ^=0 and £=2 are Raman active. This as- 
sumes that the nanoparticlc is perfectly spherical in 
shape and spherically symmetric in all of its properties. 
The basic nature of the £ = modes is much more clear 
because of their simplicity and symmetry. Consequently, 
the modes (SPH,£ = 2,n) are of primary interest when 
trying to understand Raman intensities. 

From the value of L2 ~ 0.14 in Tab. 12] the displace- 
ment of (SPH, 2,0) is mostly due to the ut term and 
not the ul term. Rs squared displacement due to the 
ul term alone is just 14% of the total. The ut term has 
zero divergence. Therefore, (SPH, 2,0) doesn't have much 
divergence. So the effect of changing density on the di- 
electric constant through the deformation potential may 
not be significant to the overall Raman intensity. 

On the other hand, based on its URS of ~ 0.66 and 



TABLE II: Other features of SPH modes for a material with 
vt/vl = 0.5. 
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UTS of ~ 0.20, the surface displacement of (SPH,2,0) 
is strongly along r and only weakly along as Fig. 
illustrates. Note that, r surface displacement changes 
the nanoparticlc shape, while 6 displacement docs not. 

The (SPH,2,1) mode differs from (SPH,2,0) in several 
ways. From the LI value of 0.2281 in Fig. 21 we can 
see that the frequency of (SPH, 2,1) depends more on vl- 
Also, L2 ~ 0.416 in Fig.Hshows that (SPH,2,1) has more 
of a ul component, even if it is still weaker than the ut 
part. But this means that (SPH, 2,1) can have much more 
divergence than (SPH, 2,0). So the deformation potential 
mechanism can modulate the dielectric constant. But it 
is very interesting to notice from the URS value of ~ 0.00 
in Tab. HH (more precisely, 0.0003) that (SPH,2,1) causes 
negligible radial movement of the surface. So (SPH, 2,1) 
barely changes the shape of the nanoparticle, as Fig. 
shows. 

(SPH, 2, 3) has strong vl dependence (LI ~ 0.8475) in 
Fig. 0] as well as a strong ul component (L2 ~ 0.766). 
So it is clear that it is SPH^. R's surface displacement is 
mostly along r and not 8 from its URS value of 0.71 and 
UTS ~ 0.04. So (SPH,2,3) will strongly affect the shape 
of the nanoparticle surface, as shown in Fig. ^ 

ul and ut take on simpler forms as r\ becomes larger. 
For large 77, the ul term has primarily radial displace- 
ment, while the ut term corresponds to displacement in 
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the 9 direction. 

For the lowest modes, the situation is qualitatively dif- 
ferent. Consider (SPH,2,0) with vt/vl = 0.5. Suppose 
to simplify this discussion we normalize the displacement 
field so that (u\u) = 1. Then (ui\ui) ~ 0.14. However, 
(u t \ut) ~ 1.85. So L3 for (SPH,2,0) is actually ~ -0.85, 
making it "ultra transverse" by that measure. It seems 
quite odd that the ut term alone has a magnitude much 
greater than that of the overall motion. The resolution 
of this puzzle is that ul and ut are not mutually orthog- 
onal with respect to the inner product of Eq. 0] In fact, 
(ut\ul) — -0.50. According to the usual vector relation, 
a-b = \\d\\\\b\\cos6 a b, the "angle" between ul and ut is ~ 
165 degrees for the (SPH,2,0) mode. This angle is nearly 
unchanged as vt/vl varies. Thus, it£ and ut are nearly 
antiparallel vectors in the function space of vector fields 
within the nanoparticle interior. It can be said, then, 
that the functional forms of u/, and ut are actually rel- 
atively similar. This is a bit of a surprise since one is 
curl-free while the other is divergence free. This angle 
between ul and ut rapidly approaches 90 degrees as r\ 
increases (i.e. for modes with higher n). 

As Fig. [21 shows, the starkness of the contrast between 
SPHl and SPHt modes is at its best when vt/vl is 
lower. For materials with high vt/vl such as Si and 
Ge, FSM modes tend more to be mixtures of SPHl and 
SPHt, especially at low r\. But the concept of longitu- 
dinally is quite applicable to materials such as Au and 
Ag. 



IV. HIGH FREQUENCY MODE 
CLASSIFICATION 

The reason for the dichotomy of SPH modes as SPHt 
and SPHl can be simply explained in the high frequency 
limit. Consider A^, the 2x2 determinant in Eq. El and 
its four matrix elements. Note that £/f] = vt/vl- So at 
high frequency, both 77 and £ are large. In that case, T\\ 
and T43 will be much larger than T13 and T41 because of 
their terms including factors of rj 2 . Consequently, Ag is 
very well approximated by T11T43. Since normal modes 
correspond to zeroes of A?, it is clear that there will be 
two sets of modes: those which are approximately zeroes 
of Tn and T43 respectively. The first group are SPHl 
and the second group are SPHt. 

The roots of Tn and T43 are plotted versus £ in Fig. El 
with lines with crosses for SPHl modes and lines with 
empty squares for SPHt modes. The lines with full cir- 
cles are the exact FSM mode frequencies. 

There are three kinds of situations where we don't ex- 
pect this approximation to be valid: (1) for low 77 (2) 
when it is not true that 77 3> £ and (3) where longitu- 
dinal (Tn) and transverse (T43) modes for a given £ are 
close - i.e. when the associated curves cross. Except in 
the previously mentioned places, the agreement between 
FSM and our approximation is quite good. The low r\ 
situation corresponds specifically to the similar prefac- 



tors of je for Tn and T43 not being large. It is apparent 
that Tn and T43 can only be useful as estimators of 
SPHl and SPHt mode frequencies when 77 3> I. This is 
confirmed from inspection of the lower right portion of 
Fig. El 

For large x, ji{x) ~ sin(a; — £^)/x. Therefore, for 
large £ and 77, the roots of Tn can be approximated by 
£ ~ ££ + (1 + ul)tt and the roots of T43 by 77 ~ £^ -\- nrpn 
where til > and tit > are integers. These lead to 
remarkably compact approximate expressions for SPHl 
and SPHr FSM frequencies in Hertz, respectively: 

(8) 

where d = 2R. These expressions are very suggestive of 
the formula for acoustic standing waves in a one dimen- 
sional system of length d. Table [H] shows the value of 
either til or nr for each mode. 

The behaviour observed in Fig. [21 becomes simple to 
explain. To a good approximation, SPH FSM modes are 
either SPHl or SPHt- This approximation is considered 
here to be good because it gives the right number of 
vibrational modes and it predicts their frequency with a 
reasonable accuracy. 

"Anti-crossing" is observed in Fig.EJeach time the vari- 
ation of the frequency of a SPHl mode crosses the one 
of a SPHt mode. In Fig.EJthere are two kinds of curves: 
horizontal lines for SPHt modes and descending curves 
for SPHl ones. Then, each time these curves come to- 
gether, an anti-crossing pattern appears for the FSM so- 
lutions. In Fig. El FSM frequencies 77 are plotted versus 
£ for a sphere made of a material which has vt/vl = 0.5. 
Because the SPHl and SPHt approximation curves are 
plotted, the anti-crossing patterns are clearly revealed. 
The continuation of Bessel functions to non-integer £ per- 
mits relationships among modes for different integer £ to 
be clearly seen. This is preferable to the common prac- 
tice of joining modes on such a graph with hand-drawn 
straight lines. 

V. DISCUSSION: 

Normal elastic waves in a solid have a longitudinal 
acoustic (LA) branch and two transverse acoustic (TA) 
branches. However, for FSM it seemed there are just two 
kinds: SPH and TOR. By classifying SPH modes into 
two kinds (i.e. SPHl and SPHt), there are now three 
categories of modes, as we would expect. 

We plot in Fig. the mean squared radial surface dis- 
placement (URS) at the surface of a 5 nm diameter silver 
nanoparticle for all SPH I = 2 modes. The magnitude 
of URS is in good agreement with the calculated Raman 
intensities 4 . (These calculations took into account the 
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non-linear dispersion of acoustic phonons in silver. As 
a result, the calculated vibration wavenumbers do not 
match.) As discussed before, (the SPH,2,0) mode is quite 
special even if we class it as a SPHt mode. It changes 
the surface shape and therefore contributes significantly 
to inelastic light scattering. Other harmonics contribute 
significantly only when their URS is large and this in 
turn is very well correlated to their SPH^ nature as can 
be seen in Fig. [2] 



CO 
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FIG. 5: Mean squared radial surface displacement (URS) as 
a function of wavenumber for (SPH/ = 2) modes of a 5 nm 
diameter silver nanoparticle (vt/vl = 0.464). 

Many experiments have observed peaks in Ra- 
man spectra attributed to acoustic phonon vi- 
brations of silver 8 ! 9 ! 10 ! 11 ! 12 ! 13 ! 14 siliconi^ and 
CdSttSei-,,! 17 ! 18 ! 19 ! 20 ! 21 These studies have regu- 
larly succeeded in observing the (SPH,2,0) mode and 
the (SPH,0,0) mode. A number of studies have seen 
(SPH,0,n) with n up to 4M. However, there has never 



been a clear indication of Raman scattering from 
(SPH,2,1) even though there have been determined 
efforts to see it . 

At the same time, (SPHl,2,?1£ = 0) seems like a strong 
candidate to have noticeable Raman scattering, since it 
has strong radial surface motion as well as a strong ul 
component that will give it stronger divergence in its in- 
terior. 

It should be noted that t = modes are always SPHl- 
That is why no full circles are plotted in Fig. Eon the T43 
root curves at I = 0. This has been the source of many 
erroneous calculations in the pasli 2 ^. 

It is often claime d 23 ' 24 ! 2 ^' 2 ^ that modes with n = 
are "surface modes" while modes with n > are "inner 
modes". The values of U2S in Tab. ITT1 show that this is a 
misconception. While this is true for I = and 1 , for I — 
2, 3 and 4 it can be seen that (SPH,^,1) has the strongest 
surface motion relative to all (SPH/,n). 

Although Tab.|n|shows URS to be zero for (SPH,2,1), 
the more precise value of vt/vl where URS is zero 
is 0.488. URS for (SPH,2,I) is only near zero for 
vt/vl close to 0.488. However, URS remains small for 
(SPH,2,1) for materials whose Poisson ratio is close to | 
which is true of many common materials. This contra- 
dicts a widespread misconception 27 ; 28 ! 29 that SPH FSM 
modes always have a radial displacement component at 
the surface. 
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